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We present a method for calculating large numbers of power spectra Ct and P(k) that accelerates 
CMBfast by a factor around 10 3 without appreciable loss of accuracy, then apply it to constrain 
11 cosmological parameters from current Cosmic Microwave Background (CMB) and Large Scale 
Structure (LSS) data. While the CMB alone still suffers from several degeneracies, allowing, e.g., 
closed models with strong tilt and tensor contributions, the shape of the real space power spectrum 
of galaxies from the IRAS Point Source Catalogue Redshift (PSCz) survey breaks these degeneracies 
and helps place strong constraints on most parameters. At 95% confidence, the combined CMB and 
LSS data imply a baryon density 0.020 < u>b < 0.037, dark matter density 0.10 < uidm < 0.32 with 
a neutrino fraction f v < 38%, vacuum density Qa < 0.76, curvature —0.19 < fik < 0.10, scalar 
tilt 0.86 < n s < 1.16, and reionization optical depth r < 0.44. These joint constraints are quite 
robust, changing little when we impose priors on the Hubble parameter, tilt, flatness, gravity waves 
or reionization. Adding nucleosynthesis and neutrino priors on the other hand tightens constraints 
considerably, requiring Qa > 0.49 and a red-tilt, n s < 1. 

The analysis allows a number of consistency tests to be made, all of which pass. At the 95% level, the 
flat scalar "concordance model" with £Ia = 0.62, uidm = 0.13, u>h = 0.02, /„ ~ 0, n s = 0.9, r = 0.1, 
h = 0.63 is consistent with the CMB and LSS data considered here, with big bang nucleosynthesis, 
cluster baryon fractions and cluster abundance. The inferred PSCz bias b ~ 1.2 agrees with the 
value estimated independently from redshift space distortions. The inferred cosmological constant 
value agrees with the one derived independently from SN la studies. Cosmology seems to be on the 
right track! 



I. INTRODUCTION 

The cosmic microwave background (CMB) is dramat- 
ically improving our knowledge of cosmological parame- 
ters 0||, although as can be seen in Figure @, the con- 
straints from CMB alone are weaker than is sometimes 
claimed. The figure shows that the adiabatic inflationary 
model that fits the CMB data best is still a crazy one with 
extreme scalar tilt n s = 1.4, curvature (f2tot = 1-3) and 
the COBE signal explained almost entirely with tensor 
fluctuations (gravity waves). 

Indeed, CMB data have now become so sensitive that 
the key issue in cosmological parameter determination is 
not always the accuracy with which CMB power spec- 
trum features (such as the position of the first peak) 
can be measured, but often what prior information is 
used or assumed (e.g., that there are no tensor fluctua- 
tions). A range of priors were explored in recent stud- 
ies including assumptions about reionization and 
gravity waves and constraints from nucleosynthesis, su- 
pernovae, large-scale structure and the Hubble constant. 
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1. The importance of including non-CMB information. 
The best fit to all current CMB data (solid curve) is still a crazy 
tilted model (n B = 1.4) for a closed universe (fitot = 1-3) with the 
COBE signal explained almost entirely with tensor fluctuations. 
The dashed curves show the tensor and scalar contributions. 
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In particular, two extensive multiparameter analyses 
included large-scale structure (LSS) information as quan- 
tified by a normalization and the so-called "shape pa- 
rameter" r (which slides a fixed transfer function side- 
ways) from galaxy power spectrum measurements and 
cluster abundances, extending earlier CMB+LSS work 

While much of the information in the galaxy power 
spectrum is indeed encapsulated in a horizontal and a 
vertical offset, all of it is clearly not. It should there- 
fore be possible to do still better by fitting directly to 
the LSS data, explicitly including the way in which each 
of the cosmological parameters affect this curve, just as 
is presently done for the CMB. This is the goal of the 
present paper. Now is a particularly exciting time to 
start doing this, since projects like the 2dF Survey and 
the Sloan Digital Sky Survey will soon produce dramatic 
improvements in LSS data quality. 

The LSS data used in this paper is the linear real 
space power spectrum of the IRAS Point Source Cata- 
logue Redshift Survey |l9| (PSCz) as measured by ]2(J 
and as shown in Figure H. The PSCz survey contains 
redshifts for 14,677 galaxies covering 84% of the sky to a 
usable depth of about 400 /i _1 Mpc. Although other large 
data sets such as the Las Campanas Redshift Survey Q] 
and the CfA/SSRS UZC redshift survey @ have compa- 
rable numbers of galaxies, the large volume of the PSCz, 
along with the careful attention paid by its authors to 
uniformity of selection, makes PSCz the most powerful 
publicly available probe of LSS at large, linear scales. 
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FIG. 2. Why LSS data adds information to the CMB. The wig- 
gly curve corresponds the best fit model from CMB alone that was 
shown in Figure 1, normalized on small scales (it has bias b = 0.5). 
The straighter curve shows Prior P5 from [Q — the difficulty in 
matching the largest scales illustrates that there is already more 
information in the curve than a normalization and a shape param- 
eter. Of the PSCz measurements shown, we opt to use only those 
in the fairly linear regime k < 0.3 h Mpc" 1 (black), discarding the 
rest. 

Afficionados may notice that the error bars on the 



PSCz power spectrum in Figure || appear somewhat 
larger than some other published measurements. This 
is because the measurements have been decorrelated pi| 
so that each plotted point represents an essentially inde- 
pendent piece of information. Having uncorrelated data 
points, or equivalently a full covariance matrix, is prereq- 
uisite for a reliable likelihood analysis. 

A longstanding obstacle to interpreting LSS measure- 
ments is the thorny issue of galaxy-to-mass bias. Local 
bias models predict that the bias factor should be con- 
stant at large, linear scales [^4-27|, and A-body exper- 
iments tend to confirm this notion |2§|-|3l|] . The sim- 
ple situation at linear scales contrasts with the non- 
linear regime, where the afore-referenced A-body ex- 
periments suggest that there is likely to be substan- 
tial scale-dependent bias. For this reason, we confine 
the analysis of the present paper to the linear regime, 
k < 0.3 /iMpc -1 . We return to this issue below. 

We will investigate how CMB and LSS constrain cos- 
mological parameters, both jointly and in separate ways 
that allow consistency checks to be made. In order to 
be able to study the effects of prior assumptions, this 
forces us to work in an 11-dimensional parameter space. 
To make this feasible in practice, we first need to de- 
velop and test a method for computing theoretical power 
spectra Ci and P(k) accurately and rapidly. 

The rest of this paper is organized as follows. In Sec- 
tion ||, we summarize our methods for computing CMB 
and LSS power spectra, saving the implementation de- 
tails and tests of their accuracy for Appendices A, B 
and C. We present our constraints on cosmological pa- 
rameters in Section III and discuss our conclusions in 



Section IV 



II. METHOD 

Our method is based on the one described in Q , but 
with a number of extensions and improvements as de- 
tailed in Appendices A, B and C. It consists of the fol- 
lowing steps: 

1. Compute power spectra Cg and P{k) for a grid of 
models in our 11-dimensional parameter space. 

2. Compute a likelihood for each model that quantifies 
how well it fits the data. 

3. Perform 11-dimensional interpolation and margin- 
alize to obtain constraints on individual parameters 
and parameter pairs. 

Our main improvement over |32| is in step 1. As de- 
scribed in Appendix B, we enhance the technique for 
accelerated CMB power spectrum calculation so that it 
becomes essentially as accurate as CMBfast itself, but 
about 10 3 times faster. We also add a simple but ac- 
curate technique to compute the grid of matter power 
spectra rapidly as described in Appendix C. In addition, 
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we improve the choice of parameters and gridding from 
|fj2f as detailed in Appendix A. For the reasons given 
there, we use the 11 parameters 

p = (t, fi k ,OA,w dm ,w b , f v ,n s ,n t , A s ,A t ,b). (1) 

These are the reionization optical depth r, the primor- 
dial amplitudes A s , A t and tilts n s , n t of scalar and 
tensor fluctuations, a bias parameter b defined as the 
ratio between rms galaxy fluctuations and rms matter 
fluctuations on large scales, and five parameters speci- 
fying the cosmic matter budget. The various contribu- 
tions fij to critical density are for curvature f2k, vac- 
uum energy cold dark matter f2 c dm, hot dark mat- 
ter (neutrinos) Q„ and baryons f2b- The quantities 
= h 2 Q,\y and u>dm = fc 2 fidm correspond to the physical 
densities of baryons and total (cold + hot) dark matter 



■ and f v = tt u /£ld m is the fraction of 



the dark matter that is hot. We assume that the bias b is 
constant on large scales but make no assumptions about 
its value, and therefore marginalize (minimize) over this 
parameter before quoting constraints on the other ten. 



III. RESULTS 
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FIG. 3. The CMB and LSS power spectra for the "concordance" 
model from Table 1. Animated versions of this figure, where the 
effect of changing one parameter at a time can be viewed, are avail- 
able at www.hep.upenn.edu/~max/concordance.html. The plotted 
model has bias b = 1.24, redshift distortion parameter /3 = 0.47 and 
hubble parameter h = 0.63. It provides an acceptable fit to all our 
data, with \ 2 ~ 96 for about 109 - 11 = 98 de grees of freedom. 

A. Basic results 

Our constraints on individual cosmological parame- 
ters arc listed in Table 1 for three cases. The best 



fit model is shown in Figure for case 3. Constraints 
are plotted in figures |i| and |] for cases 2 and 3. All 
tabulated and plotted bounds are 95% confidence lim- 
itsQ. The first case uses constraints from CMB alone, 
which arc still rather weak because of degeneracy prob- 
lems such as the one illustrated in Figure |l|. The sec- 
ond case combines the CMB information with the power 
spectrum measurements from PSCz, and is seen to give 
rather interesting constraints on most parameters except 
the tensor tilt nj^j. The third case adds three assump- 
tions: that the latest measurements of the baryon den- 
sity w b = 0.019 ± 0.0024 from Big Bang Nucleosynthesis 
(BBN) are correct |34|,|35| , that the la constraints on the 
Hubble parameter are h — 0.74 ± 0.08 |36|, and that the 
neutrino contribution is cosmologically negligible. Since 
the quoted BBN error bars are much smaller than our 
Wb grid spacing, we simply impose uj^ = 0.02. Also 
for simplicity, we take the errors on h to be Gaussian. 
The neutrino assumption is that there is no strong mass- 
degeneracy between the relevant neutrino families, and 
that the Super-Kamiokande atmospheric neutrino data 
therefore sets the scale of the neutrino density to be 
w v ~ xl0~ 4 — 10~ 3 J37|. We emphasize that this last as- 
sumption (that the heaviest neutrino weighs of order the 
root of the squared mass difference Am 2 ~ 0.07 eV 2 ) is 
merely motivated by Occam's razor, not by observational 
evidence — the best current limits on f u from other as- 
trophysical observations (see |3£| and references therein) 
are still compatible with /„ ~ 0.2. Rather, we have cho- 
sen to highlight the consequences of this prior since, as 
discussed below, it has interesting effects on other pa- 
rameters. 

For the first 7 parameters listed in Table 1, the 
numbers were computed from the corresponding 1- 
dimensional likelihood functions (these are plotted in Fig- 



* Bayesean 95% confidence limits are in general those that 
enclose 95% of the area. In this paper, we make the approx- 
imation that the boundary of the confidence region is that 
where the likelihood has fallen by a factor e~ 2 from its max- 
imum for 1-dimensional cases (such as the numbers in Table 
1) and by a factor e~ 6 ' 18 for 2-dimensional cases (such as fig- 
ures [l| and |^) . As shown in Appendix A of |52| , this approx- 
imation becomes exact only for the case when the likelihood 
has a multivariate Gaussian form. We make this approxima- 
tion to be consistent with the multidimensional marginaliza- 
tion algorithm employed here (and by most other authors), 
which is equivalent to the integration technique only for the 
Gaussian case. To give the reader a quantitative feeling for 
the importance of these issues, we also quote one-sided limis 
the on r and v in Table 1, since they have the most asym- 
metric distributions. For a detailed discussion of these issues, 
see §|. 

^The reason that we get no constraints on nt is that models 
with At = fit the data very well, for which varying n t of 
course has no effect. 
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ure [| and Figure [| for the second and third cases) . The 
best fit value corresponds to the peak in the likelihood 
function and the 95% limits correspond to where the like- 
lihood function drops below the dashed line at e~ 2 of the 
peak value. For the remaining parameters listed, which 
(except for b) are not fundamental parameters in our 11- 
dimensional grid, the numbers were computed as in |^] 
by calculating the likelihood-weighted means and stan- 
dard deviations over the the multidimensional parameter 
space. Here the tabulated best fit values are this mean 
and the limits are the mean ±2er. 
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FIG. 4. Constraints on individual parameters using only CMB 
and LSS information. The quoted 95% confidence limits are where 
each curve drops below the dashed line. 



Table 1 — Best fit values and 95% confidence limits on cosmo- 
logical parameters. The "concordance" case combines CMB and 
PSCz information with a BBN prior ui b = 0.02, a Hubble prior 
h = 0.74 ± 0.08 and a prior that /„ ~ 10 -3 . A dash indi- 
cates that no meaningful constraint was obtained. The redshift 
space distortion parameter is f3 = /(f2 m , Q^)/b, where / is the 
linear growth rate. Zi on is the redshift of reionization, to is the 
present age of the Universe and m v is the sum of the neu- 
trino masses. The values labeled as "best" are in all cases the 
ones maximizing the likelihood. For the numbers below the hor- 
izontal line, were the limits were computed from moments as de- 
scribed in the text, the corresponding mean valus are h = .53, 
Zion = 7, t = 15.6 (CMB alone), b = 1.26, h = .59, (3 = .63, 
Zion = 9, t = 13.3 (CMB+PSCz) and b = 1.10, h = .68, (3 = .51, 
Zion = 6, to = 13.4 ("concordance"). If the reader wishes to use 
some of these model pararameters for other purposes, the num- 
bers to use are thus those in the table, not the ones here in the 
caption. Since the distributions for r and /„ are quite asymmet- 
ric, we also quote the 1-sided 95% limits r < 0.22 (CMB only), 
r < 0.34, /„ < 0.35 (CMB+PSCz), r < 0.16 ("concordance"). 
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FIG. 5. Like the previous figure, but adding the "concordance" 
priors uj b = 0.02, h = 0.74 ± 0.08 and /„ ~ 10" 3 . 



B. Effects of priors 

Case 3 in Table 1 is but one from a selection of about 20 
different priors that we tried experimentally. The reason 
that we have chosen to highlight this one is that the nu- 
cleosynthesis constraint was the one that had the greatest 
impact on the results. Otherwise, the joint CMB+PSZc 
constraints were remarkably robust to prior assumption. 
Imposing priors such as flatness (Slk = 0), no tensors 
(r = 0), no tilt (n s — 1), no reionization (r = 0), and a 
reasonable Hubble parameter (we tried both h — 0.74 ±8 
at 65% and the weaker constraint 50 < h < 100 at 95%), 
both alone and in various combinations, has little effect. 
The fact that the best fit parameter values are not ap- 
preciably altered reflects that these priors all agree well 
with what is already borne out by the CMB+PSCz data: 
f2k ~ r ~ r ~ 0, and n s = 1. The fact that these priors 
do not shrink the error bars much on other parameters in- 
dicates that the PSCz has already broken the main CMB 
degeneracies. 

The nucleosynthesis prior has a greater influence be- 
cause it does not agree all that well with what the 
CMB+LSS data prefer. Although adding PSCz is seen 
to pull down the preferred baryon density slightly, reduc- 
ing the 95% lower limit from 0.024 (CMB only) to 0.020 
(CMB+PSCz), the preferred value of 0.028 still exceeds 
the BBN value. It is well-known that CMB likes either a 
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high baryon density or a red-tilt n < 1 because of the low 
second acoustic peak |l]-[i|,[59|,[lo| , and the CMB exclusion 
region to the lower right in Figure [?| illustrates this trade- 
off. Enforcing the BBN baryon value therefore shifts the 
preferred tilt-range away from the scale-invariant n s ~ 1 
case to n s ~ 0.9. Since lu\> is one of the few parame- 
ters affecting the relative heights of the acoustic peaks 
(together with Wdm, n s and, marginally, /„), eliminating 
the uncertainty in lo^ with the BBN prior also tightens 
the constraints on these other parameters. In particular, 
the link between uj^ and u>dm 1S illustrated in Figure ^. 

We found one additional prior that had a non- 
negligible effect: that on neutrinos. As illustrated in 
Figure inclusion of neutrinos substantially weakens the 
upper limits on the dark matter density. Since the neu- 
trino fraction f v has only a weak effect on the CMB, 
this effect clearly comes from LSS. A larger dark mat- 
ter density Wdm pushes matter-radiation equality back 
to an earlier time, shifting the corresponding turnover 
in P(k) to the right and thereby increasing the ratio of 
small-scale to large-scale power. Increasing the neutrino 
fraction counteracts this by suppressing the small-scale 
power (without affecting the CMB much), thereby weak- 
ening the upper limit on o->d m - Imposing the prior /„ = 
alone, without nucleosynthesis or Hubble priors, tightens 
the CMB+PSCz constraint w dm < 0.32 from Table 1 to 
u> dm < 0.19. 

Since the constraints on uJdm are tightened by fixing 
both tub (Figure ||) and /„ (Figure [)]), the "concordance" 
case in Table 1 gives quite tight constraints on the dark 
matter density. The h prior helps turn this Wdm con- 
straint into a measurement of fidmj and the measurement 
f2k ~ therefore gives an indirect constraint on the cos- 
mological constant via Q\ = 1 — £\ — ^dm — ^b- This is 
illustrated in Figure ^, where the concordance constraints 
close off the allowed region by placing a lower limit on 
fi\ ■ This lower limit on Q\ goes away if we drop either 
the Wb-prior or the /„-prior. 




FIG. 6. Constraints in the (!l m , f2A)-plane. The shaded regions 
are ruled out at 95% confidence by the information indicated. The 
allowed (white) region is seen to be centered around flat models, 
which fall on the dashed line. 



0.12 H 




Scalar spectral index n a 

FIG. 7. Constraints in the (n s , ^b)-plane. Note that PSCz not 
only shrinks the allowed region (white), but also pushes it slightly 
down to the left (the dashed line indicates the CMB-only bound- 
ary). 
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2. Ionization redshift 

In the approximation that the diffuse hydrogen in the 
Universe became fully ionized rather abruptly at a red- 
shift Zi on , this quantity is well approximated by (e.g., 
Peebles 1993) 

*on = 8.9 f^-J ^ 3 (2) 

as long as Zi on ^> 1. Although the constraints in Table 1 
agree well with what is predicted in recent simulations, 
more extreme models are seen to be ruled out. Earlier 
CMB constraints on Zion were studied in |}42|,|43f . 



FIG. 8. Constraints in the (c<Jdm> i^>b)-pl ane - As in the previous 
figure, adding PSCz prohibits high baryon solutions and allows 
slightly lower i^b-values than CMB alone. The dashed curve within 
the allowed (white) region show the sharper constraint obtained 
when imposing the priors for a flat, scalar scale-invariant model 
(Ck = r = 0, n s = 1). The dotted curve shows the effect of 
requiring negligible neutrino density (/„ ~ 0) in addition. 




FIG. 9. Constraints in the (tiAj m , /„)-plane. The shape of the 
allowed (white) region explains why the prior f v =0 tightens the 
upper limit on the dark matter density. The vertical line shows the 
CMB-only boundary before PSCz is added. 



3. Age of Universe 

The ability to place constraints in the (f2 m , r^-plane 
is allowing the age of the Universe to be predicted with 
improved accuracy The 95% confidence interval 

for our concordance case, 12.1-14.6 Gigayears, is consis- 
tent with direct age determinations from, e.g., globular 
clusters HQ. 

IV. DISCUSSION 

We have presented a method for rapid and accurate 
power spectrum calculation for large numbers of CMB 
models and used it to constrain jointly 11 cosmologi- 
cal parameters from current CMB and galaxy clustering 
data. Perhaps the most interesting results of this pa- 
per are the numbers themselves, listed in the CMB+LSS 
columns of Table 1, and their striking robustness to im- 
posing various priors. A superficial glance at the con- 
straint figures might suggest that little has changed since 
the first analysis of Boomerang + Maxima ||, or even 
since the pre-Boomerang analysis of |33] , since the plots 
look rather similar. However, whereas these earlier pa- 
pers obtained strong constraints only with various poorly 
justified priors such as no tensors, no tilt or no curvature, 
the joint CMB + LSS data are now powerful enough 
to speak for themselves, without needing any such prior 
props. 



C. Constraints on other parameters 

1. Hubble parameter 

Table 1 shows that the constraints on the Hubble pa- 
rameter h are quite weak for the CMB+LSS case. How- 
ever, the bound h < 0.78 (95%) for the concordance case 
is noteworthy since it is much stronger than (and hence 
independent of) the prior h = 0.74 ± 0.08 that was used. 



A. New public software 

Our new power spectrum calculation method accel- 
erates CMBfast by about a factor 10 3 . It is accu- 
rate to about 1-2% for Ci on all angular scales and to 
about 1% for P(k) on scales k < 0.15/iMpc -1 . Since 
this is roughly the intrinsic accuracy level of CMB- 
fast itself, there is no reason not to take advantage of 
this technique when constraining cosmological param- 
eters. A modified version of CMBfast incorporating 
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our fc-split method will be made publicly available at 
www.sns.ias.edu/~matiasz/CMBFAST/cmbfast.html. 



B. Caveats 

Let us now discuss assumptions and approximations 
that under ly our analysis. 

For both the CMB and LSS likelihood calculations, the 
percent level inaccuracies in our power spectrum compu- 
tation are likely to have a negligible effect. Indeed, the 
least accurate models tend to be wild ones that are in- 
consistent with the data in any case. On the LSS side, 
the dominant uncertainties are likely to be related to the 
measured P(k) instead. Specifically, our use of the mea- 
surements of the PSCz real-space power Q assumed 
both that linear perturbation theory was valid and that 
the bias was scale-independent on these scales. Let us 
now discuss both of these assumptions in turn. 

To assess the possibility that nonlinear effects at k ~ 
0.1-0.3 /iMpc" 1 had tainted our results, we repeated the 
entire analysis twice, discarding all P(fc)-measurements 
for k exceeding 0.2/iMpc -1 and 0.1 ftMpc -1 , respec- 
tively. The upper limit /„ < 0.38 on the neutrino frac- 
tion was weakened for the 0.2 case and went away com- 
pletely for the 0.1 case. Thus the upper limit on neutri- 
nos is sensitive to information at mildly nonlinear scales. 
The upper limit on n s was also weakened as the lever 
arm shortened, but only very slightly, from n s < 1.16 
at 0.3/iMpc -1 to n s < 1.18 and n s < 1.19 at 0.2 and 
0.1 /iMpc -1 , respectively. Other constraints were less af- 
fected. 

The redshift distortion study reported in sug- 
gested that, while nonlinear effects are visible in the 
galaxy-velocity power spectrum already at scales k w 
0.15 h Mpc" 1 , linear theory is probably a fair approxima- 
tion down to k « 0.3 h Mpc^ 1 . A subsequent study [47]] 
of the nonlinear power spectrum of PSCz has shown that 
the galaxy power spectrum is likely to be antibiased rela- 
tive to the matter power spectrum at translinear scales, a 
conclusion previously arrived at by |2^ , ^8|] . Such antibias 
tends to cancel the effects of nonlinearity in the matter 
power spectrum, making the galaxy power spectrum ap- 
pear similar to the linear matter power spectrum down to 
k w 0.3 h Mpc^ 1 . The fact that the relative bias between 
APM and PSCz galaxies is consistent with being con- 
stant, 6 A pm/&pscz « 1.15, for k < 0.3/iMpc" 1 gg^ 
adds further circumstantial evidence suggesting that the 
galaxy power spectra are near linear at these scales. 

We therefore conclude that including the PSCz data 
at k = 0.1-0.3 h Mpc -1 does not bias the results signifi- 
cantly, and that, aside from the constraint on neutrinos, 
the information from PSCz is not dominated by these 
last few bins. 

Although most theoretical work has suggested that 
the bias b is unlikely to vary much on linear scales, we 
must still be open to the possibility of scale-dependent 



bias masquerading as a cosmological effect. For instance, 
one might imagine that more luminous galaxies are more 
highly biased, and that they carry a greater statistical 
weight for the leftmost fc-bands since they remain in the 
magnitude-limited 0.6 Jy sample even at great distances. 
Such a luminosity bias could masquerade as a slight red- 
tilt n s < 1. However, three recent studies ]50|-|52^| all 
conclude that PSCz galaxies of different luminosity clus- 
ter similarly. 

On the CMB side, a long list of approximations in 
our treatment were discussed in [ p2[ , involving both the 
likelihood calculation Q and the marginalization. The 
dominant limitation is likely to be that we have not in- 
cluded the full window functions and slight band-band 
correlations of Boomerang and Maxima, but this is un- 
fortunately not possible until the relevant Fisher matrices 
are made public by the two experimental teams. 

Positive correlations between neighboring data points 
makes is easier (in terms of chi-squared cost) to shift the 
overall height of say the first peak up and down, thereby 
weakening the constraints on parameters that are mea- 
sured mainly from peak heights (wt>, w c d m , n s ). This 
effect is similar to that of calibration errors, which we 
did include. The fact that our best fit models tend to 
predict a lower first peak than much of the data illus- 
trates the effect of calibration errors and also shows that 
an extremely high first peak is hard to achieve given the 
other constraints. 

The extraction of constraints on individual parame- 
ters from the multidimensional likelihood function can 
be done in a number of different ways. Three ways of 
marginalizing are discussed in |32| (by integration, by 
maximization over the grid and by maximization over a 
smooth interpolating function), all of which have been 
used in the recent literature. An encouraging indica- 
tion that our results are insensitive to the method of 
marginalization comes from comparing the marginalized 
constraints from Table 1 with those obtained without any 
marginalization, from the above-mentioned computation 
of means and standard deviations by summing over the 
grid. These completely different methods give quite sim- 
ilar results for all well-constrained parameters. A typical 
example is the 2a range for Wd m for the "concordance" 
case, coming out as 0.110 < Wd m < 0.170 (marginalized) 
and 0.105 < Wd m < 0.172 (from moments). Even for 
r, the parameter with the most non-Gaussian likelihood, 
the 2cr "concordance" upper limits are similar: 0.161 and 
0.163, respectively. 

Finally, although we repeatedly referred to the "no 
prior" case for our CMB + LSS analysis, it is impor- 
tant to bear in mind that there is strictly speaking no 
such thing as no priors. Specifically, all our calculations 
assumed that the adiabatic inflationary paradigm is cor- 
rect. We also assumed that the dark energy was a cosmo- 
logical constant rather than some form of "quintessence" 
with a different equation of state. Finally, the edges of 
our parameter grid imposed a hard-wired top hat prior. 
This had a negligible effect on our results for all param- 
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eters except one, since the likelihood dropped to neg- 
ligible values well before reaching the boundary of our 
1 1-dimensional parameter space. The one exception in- 
volved Qa, since Figure |6| illustrates that negative Q\ 
cannot be excluded except when either BBN or SN la 
information is included. 



C. Comparison with other recent work 

Our results agree fairly well with the recent constraints 
from other groups @J^-^ • The analysis most comparable 
to ours is that of Jaffe et al. §. That study limits LSS- 
information to that incorporated in a shape parameter 
and a normalization parameter, and uses a smaller CMB 
data set, limiting the analysis to COBE, Boomerang and 
Maxima. The main effect of this culling is likely to en- 
ter on scales 50 < I < 200, covering the rise toward the 
first acoustic peaks, where Boomerang and Maxima are 
both sample variance limited and other experiments have 
covered a substantially larger sky area. In addition, || 
limit their parameter space to no neutrinos and no ten- 
sors {f v = r = 0), employ a different numerical marginal- 
ization scheme and include the above-mentioned propri- 
etary band correlations. The fact that our results agree 
so well despite all these technical differences is quite re- 
assuring, indicating that the data are now good enough 
to make the results insensitive to details of method. The 
stronger upper CMB+LSS limit on Wd m obtained in ||, 
which indirectly gives Sl\ > 0.5 as described above, is 
presumably due to their no-neutrino prior /„ = 0, as can 
be understood from Figure §. The fact that their CMB- 
only constraints are stronger than ours traces back to 
their no-tensor r = prior — this prior eliminates the 
n s —r — LUb degeneracy that we needed the full PSCz data 
to break, and automatically excludes models such as the 
one shown in Figure [lj 

Our results also agree well with those from the likeli- 
hood analysis of Kinney et al. || . Although this analysis 
has /„ = 0, a prior for nt and a limited r-range, it in- 
cludes a thorough treatment of tensor modes and maps 
out the (n s , r)-plane in detail. The study finds that quite 
blue-tilted models are allowed when r is large, precisely 
the effect that degrades our CMB-only constraints — see 
also liHU. 

D. Towards a refined concordance model 



1. Baryons 

Perhaps the most obvious one involves the baryon 
fraction. Although there is still some tension between 
BBN (preferring u h ~ 0.02 and CMB+LSS (preferring 
LJh ~ 0.03), an issue which will undoubtedly be clarified 
by improved data within a year^j, the most striking point 
is that the methods agree as well as they do. That one 
method involving nuclear physics when the Universe was 
minutes old and another involving plasma physics more 
than 100,000 years later give roughly consistent answers, 
despite involving completely different systematics, can 
hardly be described as anything short of a triumph for 
the Big Bang model. 

It is noteworthy that our addition of LSS information 
pulls down the baryon value slightly, so that a BBN- 
compatible value u>b = 0.02 is now within the 95% con- 
fidence interval. Part of the reason that that the CMB 
alone gave a stronger lower limit may be a reflection of 
the Bayesian likelihood procedure employed in this and 
all other recent papers on the topic: when a large space 
of high Ub-values are allowed, the relative likelihood for 
lower values drops. 

In all three cases listed in Table 1, the best fit model 
is consistent with the data in the sense of giving an ac- 
ceptable % 2 -value. The "CMB only" case gives \ 2 ~ 70.3 
for 87 degrees of freedom. The "CMB+PSCz" case gives 
X 2 ~ 88.7, which rises to 95.6 for the "concordance" case 
- all for w 109 — 11 degrees of freedom. The effective 
number of degrees of freedom might be a few larger than 
this, since some of the 11 parameters had little effect, 
but even taking this into account, all fits are good in the 
sense of giving reduced x 2 -values of order unity. 

Apart from BBN, our baryon value also agrees with the 
range 0.007 < fib < 0.041 inferred from a low-redshift in- 
ventory |67|] and the range 0.015 < uj^ < 0.03 at redshifts 
of a few from the Lya forest [|68|-^l| . The inferred baryon 
fraction LJh/^dm ~ 15% agrees well with that inferred 
from galaxy clusters |72|,f73| for reasonable h- values. 

It is difficult to contemplate the PSCz data in Figure |2| 
without wondering whether they show evidence for bary- 
onic wiggles in the matter power spectrum. Intriguingly, 
the location and amplitude of the wiggles in the PSCz 
power spectrum fit well to a flat, pure baryon model with 
fl A = .86 and f2 m = fl b = 0.14 (for h = 0.7), albeit with 
a large blue tilt, n s = 2, but agrees poorly with the 



It is well-known that different types of measurements 
can complement each other by breaking degeneracies. 
However, even more importantly, multiple data sets allow 
numerous consistency checks to be made. The present re- 
sults allow a number of such tests. 



* A number of possible theoretical explanations for the slight 
mismatch have been discussed in the recent literature [33| . 
Another possibility is clearly that the favored value from 
CMB will shift as data improve. Recent measurements of 
a high deuterium abundance in the interstellar medium make 
it unlikely that the standard BBN value will creep above 
cub = 0.025 [ p4^65| ], and a new QSO deuterium absorption 
study reproduces ~ 0.02 [ |S6| . 
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CMB. The PSCz data are also entirely consistent with a 
wiggle-free spectrum. 

2. Dark energy 

Another important cross-check involves the cosmolog- 
ical constant. Although the constraint 0.49 < ft\ < 0.74 
from Table 1 does not involve any supernova informa- 
tion, it agrees nicely with the recent accelerating uni- 
verse predictions from SN la This agreement 
is illustrated in Figure ^, which shows the SN la con- 
straints from fldj combining the data from both teams. 
As frequently pointed out, the conclusion Q m ~ 0.35 also 
agrees well with a number of other observations, e.g., 
the cluster abundance at various redshifts and 
cosmic velocity fields |8C]], although there is still some 
internal controversy in these two areas (see, e.g., ]8l[|). 




0.5 1 1.5 2 

p = f/b 

FIG. 10. The upper panel shows x 2 f° r the linear redshift dis- 
tortion parameter beta measured from the PSCz data. Specifi- 
cally, x 2 was computed using the galaxy-galaxy, galaxy-velocity 
and velocity-velocity power spectra but no CMB information, then 
minimized over all parameters except and Qj^. The fact that 
the resulting x 2 (f2k>^A) falls near a 1-dimensional curve with lit- 
tle scatter shows that it essentially only depends on one particular 
function of these two parameters, the linear growth rate /(£7 m , 
(equivalcntly the redshift distortion parameter f3 = f /b). The best 
fit curve (thick line) has a slower rise for high (3 than a parabola 
(thin line). The lower panel shows the corresponding likelihood 
C <x e~ x / 2 , and the la and 2a confidence limits are where this 
(thick) curve crosses the dashed lines. 



3. Bias 

A third cross-check is more subtle but equally strik- 
ing, involving the bias of the PSCz galaxies — we can 
measure it in two completely independent ways. One 
is by comparing the amplitude of the CMB and galaxy 
power spectra, which gives the constraints listed in Table 
1. The other way is via the linear redshift space distor- 
tion parameter f3 = f(Q m ,QA)/b, where f(fl m , ^a) ~ 
Sl°; 6 is the dimcnsionless linear growth rate [^2] p4j . 
We therefore remove the CMB data (which eliminates 
the /3-constraints from Table 1) and add in their place 
the two PSCz power spectra from |2(| that we dis- 
carded above. These are the galaxy-velocity and velocity- 
velocity power spectra from the stochastic bias formal- 
ism, roughly speaking corresponding to the quadrupole 
and hexadecapole of the redshift space distortions. The 
redshift distortion parameter for the PSCz galaxies was 
measured to be [3 — 0.41^ "jj |2C| ], a result in good agree- 
ment with the PSCz team's own most recent measure- 
ment, (3 = 0.39 ± 0.12 |$5|. Both of these measurements 
involved a limited marginalization over the power spec- 
trum. Here we marginalize the full likelihood function 
over all our cosmological parameters except /3. As can 
be seen in Figure [l^, this gives the la measurement 
(3 = 0.45l'^2- The fact that this agrees so well with 
the corresponding la measurement f3 = 0.51 ±0.08 from 
Table 1 for the concordance case means that a highly 
non-trivial consistency test has been passed. 

E. Concordance 

In conclusion, the simple "concordance" model in the 
last columns of Table 1 (plotted in Figure |3|) is at 
least marginally consistent with all basic cosmological 
constraints, including CMB, PSCz and nucleosynthesis. 
Specifically, as discussed above, our calculations show 
that it has passed three non-trivial consistency tests. 
Moreover our concordance model is encouragingly robust 
towards imposing a score of prior constraints in various 
combinations. Cosmology seems to be on the right track! 

The authors wish to thank Angelica de Oliveira-Costa, 
Daniel Fisher, Brad Gibson, Wayne Hu, William Kinney, 
Arthur Kosowsky, and Nikhil Padmanabhan for useful 
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NAG5-7128 and NAG5-9194, the University of Pennsyl- 
vania Research Foundation, and Hubble Fellowship HF- 
01116.01-98A from STScI, operated by AURA, Inc. un- 
der NASA contract NAS5-26555. 
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APPENDIX A: BASIC METHOD 
IMPROVEMENTS 

As mentioned, our method is based on the one de- 
scribed in [3^| , but with a number of extensions and im- 
provements as detailed below. It consists of the following 
steps: 

• Compute power spectra Ct and P{k) for a grid of 
models in our 1 1-dimensional parameter space. 

• Compute a likelihood for each model that quantifies 
how well it fits the data. 

• Perform 1 1-dimensional interpolation and margin- 
alize to obtain constraints on individual parameters 
and parameter pairs. 

Our main improvement over p2| is in step 1. We im- 
prove the technique for accelerated CMB power spectrum 
calculation so that it becomes essentially as accurate as 
CMBfast itself, simply about 10 3 times faster. We also 
add a simple but accurate technique to compute the grid 
of matter power spectra rapidly. 

1. Improved choice of parameters 

Reference p3j explored the 10-dimensional parameter 
space involving the reionization optical depth r, the pri- 
mordial amplitudes A s , A t and tilts n s , n t of scalar and 
tensor fluctuations, and various contributions Oj to crit- 
ical density. The Qi included were for curvature f2k, vac- 
uum energy Q A , cold dark matter O c( j m , hot dark matter 
(neutrinos) Vt v and baryons f2b- Since it is computa- 
tionally advantageous to work with parameters that are 
closely linked to the most important physical processes 
involved, reference |3^ | used the parameter vector 

p = (t, tt k , Q A , tu cdm , ujb, u v , n s ,n t ,A s ,A t ), (Al) 

where the physical densities u>i = h 2 £li. 

Unless il^ 1, the neutrinos left over from the early 
Universe were heavy enough to be fairly non-relativistic 
during the processes that created the acoustic peaks, and 
thereby had almost the same effect as cold dark matter on 
the CMB. The CMB power spectrum therefore depends 
mainly on the total nonbaryonic (cold+hot) dark matter 
density Wdm = w c< i m + lo v and only quite weakly on the 
hot fraction /„ = Lo v /ij dm . We therefore replace our old 
parameters (w c dm, w„) by (w,jm,/i/)- This allows us to 
accurately compute the weak /^-dependence of the scalar 
CMB power spectrum by using a coarse grid f„ =0.0, 0.3, 
1.0 and interpolating. Moreover, the tensor fluctuations 
are essentially independent of f v , so we only compute 
them for f u = 0. 

In this paper, we need to add one more parameter, re- 
lating the theoretically predicted power spectrum of mat- 
ter P(k) to that of PSCz galaxies P g (k) on large scales. 



This parameter is the bias b = \P a (k) / P(k)] 1 ^ 2 from the 
stochastic bias formalism |86|-|88f|. Although it can in 
principle depend on scale, we will assume that it is con- 
stant on the large scales that we consider [^-p^T| . We will 
make no assumptions about the value of b, however, and 
therefore marginalize over this parameter before quoting 
constraints on the other ten. In summary, we use the 
parameter vector 

p = (t, tt kl Q A ,u) dm ,u) hl f u ,n s ,n t , A Sl A tl b). (A2) 

Note that the Hubble constant is not a twelfth indepen- 
dent parameter, since 

* = V i-n k -n A - (A3) 

We wish to probe a large enough region of parameter 
space to cover even quite unconventional models. This 
way, constraints from non-CMB observations can be op- 
tionally included by explicitly multiplying the likelihood 
function C(p) by a Bayesian prior rather than being hard- 
wired in from the outset. To avoid dealing with pro- 
hibitively many models, we use a roughly logarithmic 
grid spacing for uj m , ujh and u> v , a linear grid spacing 
for f^k and f^A, a hybrid for r, /„, n s and nt, and (as 
described below) a continuous grid for A s , A t and b. 

The recent progress in CMB accuracy has been so dra- 
matic that the grids used in some recent papers ||,|| are 
already almost too sparse to accurately sample the small 
allowed regions of parameter space. We therefore modify 
the grid from p| to zoom in on the favored parameter 
ranges while still retaining some outlier points to be on 
the safe side. We let the parameters take on the following 
values: 

• r = 0,0.05,0.1,0.2,0.3,0.5,0.8 

• Ov = 0,0.1, ....,1.0 

• Q k such that Q m = 1 - O k - fl A = 0.1,0.2, 1.0 

• uj dm = .02, .05, .08, .13, .16, .20, .28, .40, .80 

• uj h = .003, .013, .016, .020, .024, .03, .04, .05, .08, .13 

• f v = 0, 0.05, 0.1, 0.2, 0.3, 0.4, 0.6, 0.8, 1.0 

• n s = 0.5, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2, 1.4, 1.7 

• m = -1.00, -0.70, -0.40, -0.20, -0.10, 

• A s is not discretized 

• A t is not discretized 

• b is not discretized 

Note that the extent of the fik-grid depends on £l A , giving 
a total of lOx 11 = 110 points in the (f£ m , iXO-plane. Our 
discrete grid thus contains 7x110x9x10x9x9x6 = 
33,679,800 models. As in p2] |, the main limitation on 
this grid size is disk space rather than CPU time. 
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2. The three basic spectra and their normalization 



It is convenient to write the two power spectra Ci and 



Pg(k) that we can measure as 



^scalar _j_ ^tensor 



a 

P g {k) = A s b 2 P(k) 



(A4) 
(A5) 



where the three basic power spectra C| calar , C| ensor 
and P{k) are all normalized consistently, correspond- 
ing to a fixed amplitude of the gravitational poten- 
tial %l> when each mode is outside horizon. The de- 
fault output of CMBfast first normalizes the CMB 
output to COBE, but the new version allows out- 
put of the raw unnormalized spectra that we need, 
^scalar depends on ( T; fik; ^ Wdm) ^ /„, „,), C~ de- 
pends on (t, J7k, ^A, w,j m , tJb, n t ) and P{k) depends on 
(f2k, W(j m , fvi n s), so we need to compute three 
separate grids of models of dimensionality 7, 6 and 6, re- 
spectively. We describe a fast and accurate way of doing 
this in Appendices B and C. 



3. Likelihoods and marginalization 

We compute the CMB likelihood exactly as in (32), i.e., 
using the first results from Maxima and Boomerang as 
well as all prior experiments |89f| (shown in Figure [l]) and 
taking into account the effect of calibration errors. 

For the PSCz galaxy power spectrum, we use only the 
band power measurements plotted in black in Figure |^, 
omitting the ones further to the right. This is a sub- 
set of the measurements from [Q that includes infor- 
mation only at scales k < CI.3ft.Mpc~ 1 to ensure that 
we stay clear of nonlinear effects. As a further precau- 
tion, we examine how the result s cha nge when this cut 
is further sharpened in Section [V B . We approximate 
the corresponding likelihood function by the multivari- 
ate Gaussian C\ ss ct e~? x , where \ is computed using 
the measurements in Figure ||. Each point in Figure || 
represents an uncorrelated measurement of the power in 
a well-defined band whose FWHM is indicated by the 
horizontal bar. In computing the likelihood, we take into 
account the detailed form of each band-power window. 

The joint likelihood function is obtained by multiplying 
the CMB and LSS likelihoods. Throughout this paper we 



marginalize over the amplitudes A s , At , and the bias fac- 



tor b. Equations (A4) and (A5) show that we can equiv- 
alently marginalize over these parameters separately for 
the CMB and LSS likelihoods before multiplying them 
together, which simplifies the calculations in practice. 



APPENDIX B: METHOD FOR COMPUTING C t 



We compute C| ensor as described in |32|: by running 
CMBfast merely for a coarser grid in Wdm and Uh (on 



which the dependence is weak) and interpolating onto 
our full grid using a regularized multidimensional spline. 
Since Cj; onsor only contributes to the first few hundred 
multipoles, it is much faster to compute than C| calar , 
which is the real challenge. 



1. The fc-space split 

The idea introduced in |5f| and |$2j was roughly speak- 
ing to compute the £ < 100 and I > 100 parts of C] calar 
separately and splice them together afterwards. The for- 
mer can be computed just as fast as C| ensor , since its 
dependence on Wdm and uj^ is weak and it is essentially 
independent of /„. The latter can be computed rapidly 
as well, since the only effect of fik and fl\ is to shift it 
sideways in a known way. To a decent approximation, the 
r-dependence can be incorporated analytically as well, as 
simply a multiplication by e~ 2r . However, since there is a 
small bump of regenerated power from the new last scat- 
tering surface (which moves to larger I as r is increased), 
we opt to include r explicitly this time — as mentioned, 
the algorithm is so fast that we are limited by disk space 
rather than CPU time anyway. In short, the high-£ part 
of (7| calar only needs to be computed on a 5-dimensional 
grid spanned by (r, w<i m , uii,, f v , n s ). Moreover, this is 
really only 4 "hard" parameters, since CMBfast treats 
multiple n s -values simultaneously with no slowdown. 

Although this approximation works well, it is typi- 
cally only accurate to 5-10% or so. The main problem 
is with the splicing itself, since projection effects alias 
power from a given physical scale to quite a broad range 
of ^-values, blurring the separation between the low and 
high grids. The new method that we present here by- 
passes this problem by making the split directly in k- 
space, where the actual physics takes place. Specifically, 
we modify CMBfast to save only the contribution to Ci 
from below or above a certain wave number fc* when it 
integrates the Boltzmann equation. 

The "low" (k < k*) contribution corresponds to fluc- 
tuations on scales outside the horizon at recombination. 
This makes it almost independent of the causal micro- 
physics that creates the familiar acoustic peaks, i.e., al- 
most independent of Wdm, and f v . Rather, it is dom- 
inated by what happens at low redshift (the late inte- 
grated Sachs-Wolfe effect, reionization, etc.). 

In contrast, the "high" (k > fc„) contribution is essen- 
tially unaffected by low redshift effects, so fik and f^A 
(or some other dark energy component that was negli- 
gible at z > 10 3 ) will merely take the pattern put in 
place at z ~ 10 3 and shift it sideways according to the 
angle-distance relationship. Suppressing the other 9 pa- 
rameters, we thus have 

Cf alar (f! k , fi A ) » w (ft k , fi A ) + C£, igh (0, 0), (Bl) 

where C = [di ss (f2k, Q\)/di ss (0, 0)]£ and di ss is the angular 
diameter distance to the last scattering surface in Mpc 
(not in /i _1 Mpc). 
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We tune the choice of differently for each model, 
choosing k* = 1.5/d s hor where d s hor = / c 8 drj is the 
sound horizon at decoupling (here c s is the sound speed 
and 77 denotes conformal time). In ^-space, this value of 
/c* corresponds approximately to the place were we did 
the splitting in |32|, i.e., to the early rise of the first 
acoustic peak (£ = 100 for flat models, higher/lower I for 
open/closed models). 




1 10 100 1000 



FIG. 11. Three examples of our method for computing 
CMB power spectra. The three panels show a fiat A-model 
(Ofc = 0, Oa = 0.7 — top), a moderately open model 
(Qk = 0-4, Qa = 0.3 — middle) and a moderately closed model 
(Qk = — 0.2.£7a = 0.9 — bottom). All three models have 
ui dm = 0.1225, u b = 0.0245 and J v = 0. The model cf gh used in 
the k > A;* calculation had Slj. = Ha = 0, and is simply shifted side- 
ways differently in each panel, whereas the (k < fc*) spectra (C^ ow ) 
were computed separately for each model. The solid line shows the 
full CMBfast calculation while the long dashed line shows the result 
of our new method, i.e., the sum of the two dashed curves. 

Figure p"l| shows an example of this splitting technique. 
We show three panels with a flat A model, a moderately 
open model and a moderately closed models. Each panel 
has four curves. The solid line is the model calculated 
fully with CMBfast from scratch. The long dashed line 
shows the same model calculated with our new technique, 
and is seen to differ by less than 2% across the spectrum. 
For completeness we also show the spectra for k < k* 
and k > fc* that were added in each panel. Note that 
the k > fc* curve is the same for the three panels, merely 
shifted sideways by different amounts to match the an- 
gular diameter distance. 

We do not make a sharp cut at k*. Rather, to avoid 



numerical problems, we use a soft cut defined by the func- 
tion 

w{k) = — It?* ( B2 ) 

Specifically, when computing C] ow and C^ lg , we mul- 
tiply the primordial k power spectrum by w{k) and 
[1 — w(k)], respectively. Note that even a sharp fc-cut 
would result in a fuzzy ^-cut, since projection effects alias 
a given A:- value onto a range of ^-values. 

2. Testing the Ci accuracy 

To test the accuracy of our method, we drew a random 
sample of ~ 10 3 of the models from our final grid and re- 
computed them from scratch with CMBfast. We also 
added about 10 2 models to the test sample by hand that 
we suspected might be particularly troublesome. The re- 
sults are shown in Figure [l2[ As can be seen, the median 
accuracy is better than 2% for most ^-values. It should 
be noted that this is a test not only of the A:-space split- 
ting technique, but of our full pipeline, which includes 
several steps of interpolation. 

One interesting thing to point out is that our median 
error is substantially lower than our mean error. This 
occurs because there are a small number of outlier mod- 
els where we find significantly larger errors. We should 
first note that our worst model is never more that 30% 
off for £ < 1000. An examination of the worst models re- 
veals our main sources of inaccuracies. Our main source 
is percent level errors is the calculation of the angular 
diameter distances, which lead to a small relative shift 
between the spectra. In models with sharp peaks, this 
can lead to large relative errors although the two curves 
are very similar to each other. Thus this is quite a benign 
error that can be further improved by a better calcula- 
tion of the angular diameter distance but that has no 
effect for the likelihood of current data which has win- 
dow functions that are much wider than this small shift. 

The second source of error is in the region were we 
combine the high and low k spectra. Usually our worst 
models (which are off by less that 15% in this region) have 
a significant ISW contribution at low I, so that this ISW 
contribution is still significant for the high k wavelengths 
but were not included because we always shift flat fi m = 
1 models. Although this source of error is inherent to our 
method, it should not be a source of concern in practice 
because models with such a large ISW contribution are 
already ruled out by the data; in particular, they are 
inconsistent with the rather flat low I power spectrum 
seen by COBE. 

Finally there were errors that could be traced to the 
coarseness of our grid. Specifically our low grid had only 
three values of r. Interpolation errors lead to differences 
(around 10% for our worst models) for some high r mod- 
els. These errors could be trivially reduced by refining 
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our r-grid although the high r models that are inaccurate 
are disfavored by the data. 
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We compute G numerically using the publicly avail- 
able GrowA package Q. There are excellent packages 
of fitting formulae available for rapid computation of the 
transfer function T |)(]|)l|], but unfortunately none of 
these currently handle the general case that we need. 
Since T(k) depends on merely three parameters, CPU 
time is not an issue and we simply compute it numeri- 
cally using CMBfast §§ . 



4. Testing the P(k) accuracy 



FIG. 12. Median differences between the results of directly com- 
puting a model with CMBfast and the results in our grid as a func- 
tion of I. The top panel shows absolute differences while the lower 
panel shows relative errors. 

APPENDIX C: METHOD FOR COMPUTING P(jfc) 



3. The approximation 
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When normalized according to equation ( |A5| ) , the mat- 
ter power spectrum P(k) depends on the six parameters 
(fik, Qa> o^dmi w b) fv, n s ). However, this dependence can 
be approximately factored as 



P(k; n k , fi A , uj dm , u h , f v , n s ) 

G(n k ,n A ) 2 T(hk;uj dnl ,u; h , f„) 2 ( f- 

fc 



(C3) 



where h is given by equation (A3). Here G is the growth 
factor from linear perturbation theory |32] |84) and T is 
the transfer function normalized so that T = 1 for k = 0. 
We use ko = 0.05 Mpc -1 to match the tilt convention of 
CMBfast. It is important to note that the wave num- 
ber entering in T is measured using physical distance 
units Mpc -1 whereas that entering in P is measured in 
astronomical distance u nits (/iMpc -1 ). The approxima- 
tion given by equation ( |C3| ) becomes exact for the case 
f v = 0, and we will quantify its accuracy for the neutrino 
case below. 



FIG. 13. Differences between the results of directly computing 
a model with CMBfast and the results in our grid as a function of 
wavelength. The error distribution of each k is illustrated by its 
median and 90th percentile. The top panel shows absolute differ- 
ences while the lower panel shows relative errors. 

To test the accuracy of our method, we drew a random 
sample of 10 3 of the models from our final grid that were 
not ruled out at more than 5a and recomputed them 
from scratch with CMBfast. The results are shown in 
Figure |l3|. We see that the median accuracy of our P(k)- 
method is better than 1% for k < 0.15 and never gets 
worse than about 1.4% over our range of interest. We also 
tested 10 3 random models from the full parameter space 
(without the 5cr cut on unphysical models), obtaining 
median errors similar to the 90% cur ve in Figure [l3| 

The only reason that equation (C3) is not exact is that 
neutrinos affect the growth rate of fluctuations at late 
times when Ok and Oa become important. f v there- 
fore cannot be separated completely from G. On the 
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other hand, /„ cannot be absorbed into G either, since 
it suppresses only small-scale fl uctu ations. Fortunately, 
Figure [l|| shows that equation ( p3| ) is nonetheless very 
accurate in practice, breaking down badly only on scales 
smaller than are relevant to our present analysis. This 
is because at the low redshifts where J7k and f^A be- 
come important, the neutrino free-streaming scale below 
which fluctuations are suppressed is below a few /i _1 Mpc 
if m„ > leV. 

Even for f v = when equation (|C3| ) is strictly speaking 
exact, the practical implementation can cause small in- 
accuracies. We found our 0.6% "noise floor" seen in Fig- 
ure |l3| to be due to the horizontal shifting of the transfer 
function given by h, which we accomplished with a cubic 
spline. This effect is also responsible for part of the rise 
in relative errors towards small scales, where high baryon 
models have pronounced wiggles. It should be possible to 
eliminate this problem by computing the input transfer 
function grid with more finely spaced fc-values, able to 
oversample all baryonic wiggles. 
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